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Vortex Fluctuations in the Critical Casimir Effect of Superfluid and Superconducting Films 
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Vortex-loop renormalization techniques are used to calculate the magnitude of the critical Casimir forces in 
superfluid films. The force is found to become appreciable when size of the thermal vortex loops is comparable 
to the film thickness, and the results for T < Tc are found to match very well with perturbative renormalization 
theories that have only been carried out for T > Tc. When applied to a high-Tc superconducting film connected 
to a bulk sample, the Casimir force causes a voltage difference to appear between the film and bulk, and estimates 
show that this may be readily measurable. 
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Critical fluctuations in a finite-size superfluid lead to a free- 
energy difference between the finite-size system and the bulk. 
If there is a connection between the two, such as between 
a saturated helium film and a bulk liquid reservoir, forces 
(known as Casimir forces, in analogy with size-limited elec- 
tromagnetic fluctuations) will develop that lead to a thinning 
of the film in the vicinity of the superfluid transition, an ef- 
fect which has been observed experimentally [1, 2]. Existing 
theories of the effect are incomplete: perturbative e-expansion 
theories [3] are only able to calculate the force in the non- 
superfluid region T > Tc, and the maximum amplitude pre- 
dicted for Dirichlet boundary conditions is about a factor of 
50 times smaller than the observed maximum [2]. The super- 
fluid regime T < Tc is apparently far more difficult for the 
perturbation theories, and such calculations have not yet been 
attempted. 

We show here that vortex excitations [4, 5] are the source 
of the critical fluctuations giving rise to the critical Casimir 
force, and that vortex renormalization techniques provide a 
very simple means of calculating the force in the superfluid 
phase T < Tc- The force becomes appreciable when the 
size of the thermally-excited vortex loops become compara- 
ble to the film thickness, and the results for periodic boundary 
conditions match very well with the perturbation theories at 
Tc- When the loops become larger than the film thickness 
there is then a crossover to two-dimensional (2D) Kosterlitz- 
Thouless (KT) vortex pairs, and this leads to a prediction that 
the KT superfluid transition will take place at a temperature 
only slightly higher than the point where the Casimir force 
begins to be measurable: in helium films nearly the entire dip 
in thickness will occur in the normal state above Tkt- 

We also propose that an analogous Casimir force should ap- 
pear at the junction between a high-Tc superconducting film 
and the bulk superconductor. In this case the force will take 
the form of an electrical potential difference appearing be- 
tween the film and bulk, due to a transfer of Cooper pairs 
from the film to the bulk that balances the Casimir energy 
difference. Rough estimates show that this may be a readily 
measurable voltage (microvolts), and could be a useful probe 
of the high-Tc superfluid transition. 

The difference in free energy per unit area between the film 
of thickness i and the bulk is given by 6F ~ L{fi,—ff) where 
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FIG. 1 : Free-energy difference between bulk and film as a function 
of the film thickness L, for several reduced temperatures t — {Tc — 

T)/Tc. 



fh and ff are the free energies per unit volume of the bulk and 
film. In the vortex-loop renormalization scheme these free 
energies [5] can be written as an integral over the average loop 
diameter a, and the difference is then 



SF 



^4 r f-Vexp(-C/(a)/fcsr)- 
flo Ji3L \ao/ ao 



(1) 



where ao is the bare core diameter ( = 2.53 A for helium pa- 
rameters [6]), U{a) is the renormalized loop energy [4], and 
j3L is the maximum loop size in the film; a comparison [6] 
with a finite-size path-integral Monte Carlo simulation gave /3 
= 0.75. In Ref. [4] it was noted that from the scaling relation 
for the superfluid density the Boltzmann factor can be written 
in the form 



exp(-t/(a)/fcBT) 



3ao 

47^3 



d_ 

da 



(2) 



where Kr = h'^ Psao/'rn^kBT is the dimensionless superfluid 
density. Figure 1 shows an evaluation of Eq. 1 as a func- 
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FIG. 2: Free energy scaling function versus L/^, where ^ is the bulk 
correlation length, taken to be positive for T < Tc and negative for T 
>Tc. 



tion of the film thickness L, using the loop recursion rela- 
tions in Ref. [4], for several reduced temperatures near Tg. 
Very close to Tc there is a crossover from exponential to al- 
gebraic decay in L, since at Tc the asymptotically exact so- 
lution of the recursion relation is Kr = Do [ao/a) where 
Do = 0.3875 is a universal constant [4, 5]. Inserting this 
into Eqs. 1 and 2 gives SF/ksT = A/L^ precisely at Tc, 
where A = — l/(47r^£'o/3^)- This is exactly the form ini- 
tially predicted by Fisher and DeGennes [7] from scaling ar- 
guments, where the universal constant A is known as the 
Casimir amplitude. With /3 = 0.75 we find A = -0.155, in 
very reasonable agreement with the value A = -0.20 found 
in the e-expansion results for periodic boundary conditions 
[3]. As in Ref. [6] our calculation is equivalent to periodic 
boundary conditions, since we assume the superfluid density 
is a constant across the film, with no variation at the wall or 
free surface. Figure 2 shows the free energy scaling function 
= L'^ SF/ksT, plotted as a function of the scaling variable 
L/^ where ^ = ao/Kr is the bulk correlation length, the size 
of the largest loops being thermally excited. The solid line 
is the vortex-loop result from evaluating Eq. 1, and where for 
plotting purposes we have taken ^ positive for T < Tc- When 
the maximum loop size becomes comparable to the film thick- 
ness the free energy difference decreases rapidly, with a finite 
slope at Tc where L/^ = 0. The dotted curve is the e-expansion 
result [3] for periodic boundary conditions, normalized by the 
ratio of the Casimir amplitudes, 0.155/0.20. It is clear there is 
good agreement between the two calculations, with the finite 
slope at Tc matching well. 

The Casimir force Kc leading to the film thinning is the 
derivative of the free energy difference, and is most conve- 
niently written [3] in terms of a scaling function ^, Kc = 
-d 5F/dL = ksT'd/L^. Figure 3 shows the results for the 
scaling function, which can be extracted from the experimen- 
tal data as in Ref. [2]. This cannot be directly compared with 



the experimental results, however, since the superfluid den- 
sity in a real film falls to zero at the boundaries, and hence 
Dirichlet boundary conditions rather than periodic should be 
appUed. The depression of the superfluid density at the sur- 
faces also has the effect of shifting the transition temperature 
downward from the bulk critical temperature T\ by an amount 
dependent on L, Tc = Tc{L). It should be possible in further 
work with the vortex-loop theory to account for the depressed 
superfluid density by calculating the excess loop density near 
awaU[8]. 

A further effect which must be taken into account is the 
crossover from 3D to the 2D KT transition. In the above cal- 
culation the iterations in the film have been stopped when the 
loop size reaches the film thickness, but actually at that point 
the loops intercept the boundaries and turn into vortex pairs 
at longer length scales. Adding these excitations to the above 
Casimir calculation is easily carried out [8] by matching the 
loop recursion relations to the Kosterlitz recursion relations[9] 
for the vortex pairs at the crossover length PL. The areal su- 
perfluid density rr,, is related to ps resulting from the loop re- 
cursion relations by as = PsL, and the fugacity of the pairs is 
set equal to the loop fugacity at that length scale. The core 
diameter of the pairs is Oc, the effective core diameter of the 
loops at the crossover[4, 5], which is proportional to the cor- 
relation length, and hence is of the order of the film thickness 
L. 

With these inputs the KT recursion relations are then it- 
erated to macroscopic length scales. The superfluid density 
jumps to zero at the temperature Tkt indicated in Fig. 3. Tkt 
is a function of the scahng variable L/^, as shown previously 
in Ref. [8] where agreement was found with finite-thickness 
scaling of the KT transition[10]. Since the film thinning also 
begins to occur when the correlation length is comparable to 
the film thickness, it is not surprising that the KT transition 
occurs close to the onset of the film thiiming, with nearly all 
of the film thinning occurring in the normal state above Tkt- 
By also iterating the Kosterlitz free energy along with the re- 
cusion relations, the 2D contribution to the Casimir forces can 
be evaluated. We find that for thick films (L/oq ~ 100) this is 
negligible compared to the contribution from the loops, since 
the Kosterlitz free energy for the vortex pairs is proportional 
to aj^ w reducing the pair free energy per unit area by a 
factor of (flo/L)^ compared to that of the loops. An important 
experimental check of these predictions for the KT onset point 
will be to measure the superfluid density simultaneously with 
the film thinning and check whether Tkt is coincident with 
the onset of thinning. Although again the present results apply 
only for the case of periodic boundary conditions, it is likely 
that both Tc and Tkt scale in just the same way for Dirichlet 
boundary conditions. 

A similar Casimir force should also be present near the 
phase transition of a high-Tc superconducting thin film con- 
nected to a bulk sample of the same material. There is now 
considerable theoretical and experimental evidence that the 
high-Tc transition is also a vortex-loop transition entirely sim- 
ilar to that of helium [4, 1 1, 12]. To give a concrete example of 
this we show in Figure 4 a fit of the vortex-fluctuation theory 
of Shenoy and Chattopadhyay [11] to recent experimental data 
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FIG. 3: Casimir force scaling function i? versus L/^. 



[13] for the superfluid fraction of a BiaSraCaCuOz (BSCCO) 
epitaxial film of thickness 610 A. This film (labeled curve B 
in Ref. [13]) is sUghtly overdoped, with a of 84.9 K. At 
low temperatures vortices are not thermally excited, and the 
only excitations that affect the superfluid density are the nodal 
quasiparticles. Fitting to the data in the range from 5 to 20 K 
gives a "bare" superfluid fraction p°/p= 1 — AT"^ with A 
= 7.5 X 10~^ K~^, and this provides a good description up to 
about 60 K where the vortices take over. The decrease of 
the superfluid fraction is the form expected for d-wave quasi- 
particles in the dirty Umit [14]. 

Above 60 K the initial vortex excitations are single-layer 
pancake-antipancake pairs[l 1]. These are particularly impor- 
tant in BSCCO due its strong anisotropy, with an anisotropy 
factor at small length scales 70 = C[i/^-L = 50 [15], where 
^11 (~ 25 A) and are the correlation lengths in the direc- 
tions parallel and perpendicular to the CuO planes. Modified 
KT recursion relations that include the linear potential from 
the Josephson cores are iterated from an initial pair separation 
ao = £,\\ as in Ref. [1 1], except that the bare superfluid fraction 
is taken to be that of the nodal quasiparticles above. The only 
fitting parameter is the core energy of a pancake vortex; the 
fit to the data gives a value of the core energy over the bare 
superfluid density I Ko = l.5,m units ksT. This is quite 
comparable to the value of 2.2 found in helium films [1 6] . 

At the length scale Tq = ^oao there is then a crossover to 
vortex loops, which initially are quite elhptical due to the 
strong anisotropy, as shown in Fig. 4. The starting fugacity 
and superfluid density of the loops are matched to those of the 
2D pairs at and for a bulk sample the 3D loop recursion 
relations are then iterated to scales greater than the coherence 
length, where the superfluid density becomes constant. The 
result for a bulk sample is shown as the solid line in Fig. 4, in 
excellent agreement with the film data at least to within a few 
degrees of T^- Beginning about 2 K from Tc the recursion re- 
lation for the anisotropy factor[l 1] shows that it scales rapidly 
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FIG. 4: Superfluid fraction in a BSCCO film, with data from 
Ref. [13]. The dashed curve is the background from the nodal 
quasiparticles, the solid curve that from single-layer vortex pairs and 
loops, and the dash-dot line indicates the location of the KT jump. 
The schematic figure above the plot indicates the progression with 
temperature of the 2D single-layer pairs, 3D anisotropic loops, and 
the crossover back to 2D vortex pairs. 
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FIG. 5: Schematic of an experiment to measure the voltage differ- 
ence between film and bulk, using two Josephson junctions. 



towards one, so that purely circular loops are excited at large 
length scales. In this limit where psj p < 0.1 the superfluid 
exponent becomes[5] v = 0.6717, the XY model value. The 
temperature range where XY critical behavior can be observed 
in BSCCO is considerably smaller than in other high-Tg ma- 
terials because of its large anisotropy factor. 

The reason the film data falls away from the predicted line 
for the bulk is hkely due to a crossover to the 2D KT transi- 
tion, which will occur when the loop size becomes compara- 
ble to the film thickness, as discussed above for helium films. 
For this case the superfluid density and fugacity are matched 
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at the scale 0.5 L, to better account for the depressed super- 
fluid density at the surfaces in the experiments. The 2D KT re- 
cursion relations are then iterated to larger length scales. The 
KT transition occurs quickly after the crossover; the dash-dot 
line in Fig. 4 shows the universal jump of the superfluid den- 
sity to zero, which coincides well with the deviation of the 
film data from the bulk curve. The experimental data is not 
expected to show a sharp jump to zero, since it is taken at a 
rather high frequency of 80 kHz. It is well known that this 
leads to a finite-frequency broadening [10, 17] quite similar 
to the behavior seen in the data here. 

The critical Casimir free energy difference between a 
BSCCO film and bulk can now be calculated in just the same 
manner as for helium films. The quasiparticle and single-layer 
pair contributions are the same in the film and bulk and cancel, 
and the KT pairs give a negligible contribution. The difference 
is entirely from the loop size being cut off by the film thick- 
ness, and the result for periodic boundary conditions is essen- 
tially the same as Fig. 2. The free energy difference leads to a 
transfer from the film to the bulk of 5N Cooper pairs per unit 
area of the film, and this gives rise to an electrical potential 
difference 5V that appears at the interface between the film 
and bulk. Since the total electrochemical potential ji + eV 
must be a constant throughout the system, the free energy 
difference is proportional to the electrical potential difference 
5V = Vb- Vf, 

6F = - Hf) SN = -e 6V 5N (3) 

But the number of Cooper pairs that transfer must also be pro- 
portional to the voltage difference[19], 

6N = -2N{Ef) L{ij,b- Hf) =2N{Ef) LeSV (4) 

where N(Eir) is the density of single-spin electronic states at 
the Fermi surface. Combining these gives a prediction for the 
voltage difference dV = y/-6F/2e^N{EF) L. This can be 
crudely estimated for a BSCCO film 600 A thick by taking a 
maximum SF/ksTc w A/L^ with A w -0.15 and « 85 



K, N(Ef) « 1 X 10^*^/ Jm3, yielding a maximum 5V « 30 /iV. 
This is the estimate for periodic boundary conditions and not 
Dirichlet, but we note that the experimental results in helium 
films appear to give an even larger value of A than used here. 

It may be possible to measure this voltage with the 
schematic experiment shown in Fig. 5, where two high-Tc 
Josephson junction contacts equilibrate with the chemical po- 
tential in the film and bulk, and the voltage difference SV 
across them can be measured with a low-Tc SQUID voltmeter. 
This technique has been used to measure voltage differences 
in low-Tc superconductors [18] to a resolution of better than 
nanovolts, and it should be possible to apply the same tech- 
nique to the Casimir voltage measurement. The electrodes of 
the two Josephson junctions need to have a higher Tc than the 
bulk or film; they can either be a different high-Tc material, or 
the same material but closer to optimal doping. Observation 
of this voltage would give a sensitive new probe of the high- 
Tc phase transition, and would allow measurements of the per- 
pendicular coherence length near the transition, since the volt- 
age is only appreciable when the coherence length equals the 
film thickness. In summary, we have shown that the critical 
Casinoir effect in superfluid and superconducting films results 
from the finite-size Umitation of vortex loop fluctuations in the 
film. This is the first calculation of the effect in the superfluid 
phase; the ease of evaluating Eq. 1 using the loop renormal- 
ization should be contrasted with the extreme complexity of 
traditional perturbative methods that prevent the calculations 
from being carried out below Tc. For periodic boundary con- 
ditions the loop result matches with the perturbation theories 
at Tc, but a full comparison with experiments will require ex- 
tension of the theory to Dirichlet boundary conditions. 
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